GRB 060313: A New Paradigm for Short-Hard Bursts? 
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ABSTRACT 



We report the simultaneous observations of the prompt emission in the 7- 
ray and hard X-ray bands by the Swift-BAT and the KONUS-Wind instruments 
of the short-hard burst, GRB 060313. The observations reveal multiple peaks 
in both the 7-ray and hard X-ray bands suggesting a highly variable outflow 
from the central explosion. We also describe the early-time observations of the 
X-ray and UV/Optical afterglows by the Swift XRT and UVOT instruments. 
The combination of the X-ray and UV/Optical observations provide the most 
comprehensive lightcurves to date of a short-hard burst at such an early epoch. 
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The afterglows exhibit complex structure with different decay indices and flar- 
ing. This behavior can be explained by the combination of a structured jet, 
radiative loss of energy, and decreasing microphysics parameters occurring in a 
circum-burst medium with densities varying by a factor of approximately two on 
a length scale of 10 17 cm. These density variations are normally associated with 
the environment of a massive star and inhomogeneities in its windy medium. 
However, the mean density of the observed medium (n ~ 1CT 4 cm 3 ) is much less 
than that expected for a massive star. Although the collapse of a massive star as 
the origin of GRB 060313 is unlikely, the merger of a compact binary also poses 
problems for explaining the behavior of this burst. Two possible suggestions for 
explaining this scenario are: some short bursts may arise from a mechanism that 
does not invoke the conventional compact binary model, or soft late-time central 
engine activity is producing UV/optical but no X-ray flaring. 

Subject headings: gamma-rays: bursts 

1. Introduction 

Gamma-ray bursts (GRBs) are generally classified into one of two categories: long-soft 
and short-hard (Kouveliotou et al. 1993) 1 . Short GRBs are bursts with T 90 < 2 s (where 
T90 is the time in which the cumulative counts above background are between 5% and 95% 
of the total) and a hardness ratio that is typically larger than the one of long GRBs. Short 
bursts also exhibit an initial spike with insignificant spectral evolution at energies greater 
than ~ 25 keV (Norris & Bonnell 2006). 25-30% of all GRBs in the Burst and Transient 
Source Experiment (BATSE) catalog (Fishman et al. 1994) are short-hard bursts (SHBs; 
Kouveliotou et al. 1993), while 12±4% of the Swift (Gehrels et al. 2004) sample are SHBs. 
It has been proposed that SHBs are the result of a compact binary merger such as a double 
neutron star (DNS; Paczyhski 1986; Eichler et al. 1989; Narayan, Paczyhski, & Piran 1992; 
Fryer, Woosley, & Hartmann 1999; Rosswog, Ramirez-Ruiz, & Davies 2003; Lee, Ramirez- 
Ruiz, & Granot 2005; Panaitescu 2006a) or a neutron star and a black hole (NS-BH; Lattimer 
& Schramm 1976; Paczyhski 1991; Narayan, Paczyhski, & Piran 1992; Mochkovitch et al. 
1993; Fryer, Woosley, & Hartmann 1999; Lee, Ramirez-Ruiz, & Granot 2005; Panaitescu 



1 Two additional classes of bursts have been proposed: an intermediate class (Horvath 1998; Mukherjce 
et al. 1998; Horvath et al. 2006) which falls between the short and long burst classes with respect to Tgo and 
has a hardness ratio softer than long bursts, and a very short class (Cline et al. 2005) which has a very hard 
spectrum and T 90 < 100 ms. 
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2006a). Since the natal kick of a neutron star is ~ 200 — 1000 kms -1 , and since the typical 
life time of compact binaries is ~ 0.1 — 1 Gyr (Panaitescu 2006a), a compact binary could 
wander anywhere from 20 kpc to 1 Mpc in the time it takes for the binary to merge. At 
these large distances from the neutron star's stellar nursery, a lower density external medium 
is expected. 

Understanding and testing these different short GRB models has been difficult due to 
the large time delay in localizing the burst afterglows. The Interplanetary Network (IPN; 
Hurley et al. 2005) has triangulated a few short GRBs on timescales of days (Bloom et al. 
2006); however, on these timescales, the afterglow has long since faded. Since the launch of 
the Swift satellite, eleven SHBs 2 have been localized in under an hour. Four of the eleven 
bursts, GRBs 050202 (Tueller et al. 2005), 050906 (Krimm et al. 2005), 050925 (Holland 
et al. 2005), and 051105A (Mineo et al. 2005), were localized in the hard X-ray band by 
the Swift Burst Alert Telescope (BAT; Barthelmy et al. 2005a) but with no corresponding 
X-ray, optical, or radio detections; three additional bursts, GRBs 050509B (Gehrels et al. 
2005; Bloom et al. 2006), 050813 (Retter et al. 2005; Morris et al. 2005), and 051210 (La 
Parola et al. 2006), were localized by the BAT with corresponding Swift X-ray Telescope 
(XRT; Burrows et al. 2005a) detections but no optical or radio detections; the remaining 
four bursts, GRBs 050724 (Barthelmy et al. 2005b; Berger et al. 2005; Campana et al. 
2006; Grupe et al. 2006; Vaughan et al. 2006), 051221A (Parsons et al. 2005; Burrows, 
Capalbi, & Grupe 2005; Grupe, Burrows, & Patel 2005; Roming et al. 2005a; Soderberg 
et al. 2006; Burrows et al. 2006), 051227 (Barbier et al. 2006), and 060313 (Pagani et al. 
2006a), were localized by the BAT and XRT with corresponding Swift Ultra- Violet /Optical 
Telescope (UVOT; Roming et al. 2005b) and/or ground-based detections. Additional rapid 
observations of GRBs 050709 (Fox et al. 2005; Villasenor et al. 2005; Covino et al. 2006) and 
060121 (Arimoto et al. 2006; Levan et al. 2006) were made with the High Energy Transient 
Explorer (HETE-2). These rapid observations have broadened our understanding of short 
GRBs and provided strong evidence that the origin of short bursts is different from that 
of long bursts and is consistent with the merger of compact objects (Gehrels et al. 2005; 
Barthelmy et al. 2005b; Fox et al. 2005; Villasenor et al. 2005). 

However, the data for GRB 060313 indicates that there are problems with the model 



2 It has been noted in the literature that there have been other well localized short bursts which have 
soft spectra (Bloom et al. 2006). Recent work (Levan et al. 2006; Sakamoto et al. 2006) has suggested that 
Swift and HETE-2 short bursts have softer spectra than the BATSE bursts. They find that the hardness 
ratio is in-between that of the BATSE short and long GRBs, but the Swift short bursts are typically harder 
than the long ones. Hereafter, we assume the Swift and HETE-2 short bursts fit into the hard classification 
scheme and refer to short bursts as those bursts with short-hard spectra. 
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for a compact binary merger, at least for this burst. Two possible explanatons for the 
observed behavior are: some short bursts may be produced by a different mechanism than 
that invoked by the standard compact binary model, or late-time central engine activity is 
generating late internal shocks from which low energy (UV/optical) flares arise. Here we 
present the broad-band nature of the short GRB 060313 and its afterglow based on Swift 
BAT, XRT, and UVOT as well as KONUS-Wind (Aptekar et al. 1995) data. This is only 
the fifth short burst with a reported optical afterglow and the first with UV detections and 
optical flaring. It is also the hardest burst in the Swift sample. In addition, the optical 
and UV detections made by the UVOT of this burst are the earliest (< 80 s) optical/UV 
detections of a short burst to date. The combination of the UVOT data with the BAT, 
KONUS, and XRT data provide the most comprehensive lightcurves to date of a SHB at 
such an early epoch and affords a unique opportunity to probe the mechanism behind this 
short GRB. 



2. Observations & Data Analysis 

On March 13, 2006, at 00 : 12 : 06.484 UT, the BAT triggered on GRB 060313 (Pagani 
et al. 2006a). The timing analysis hereafter is referenced from the BAT trigger time (T ). 
The initial BAT light curve has two slightly overlapping peaks with a total duration of 
~ 1 s. The peak count rate was approximately 90,000 counts s _1 (15 — 350 keV) occurring 
at T + 0.5 s. 

The UVOT began a finding chart exposure 78 s after the trigger. The afterglow was 
discovered during ground processing. The XRT began its autonomous sequence of observa- 
tions of the GRB field at 00 : 13 : 24 UT, 79 s after the BAT trigger. The XRT on-board 
centroiding algorithm could not converge on a source in the image, but a fading source was 
detected in the ground-processed data (Pagani et al. 2006b). The first ground based detec- 
tion was made at 01 : 28 UT with the VLT and FORS2 telescopes (Levan & Hjorth 2006). 
No radio source was detected at the VLT position (Soderberg & Frail 2006). 

Hereafter, we adopt the notation -F(z/, t) cx t~ a v~P for the afterglow flux as a function of 
time, where v is the frequency of the observed flux, t is the time post trigger, f3 is the spectral 
index which is related to the photon index r {(5 = T — 1) , and a is the temporal decay slope. 
We also adopt the convention q x = 10 x q in cgs units. A flat cosmological constant dominated 
cosmology with the parameter values H = 70 kms -1 Mpc -1 , VL M = 0.3, and f2 A = 0.7 is 
assumed. 
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2.1. BAT Analysis 

Ground analysis (Markwardt et al. 2006) of the BAT data determined that Tgo for the 
15 — 350 keV band is 0.7 ±0.1 s (estimated error includes the systematics) with a fluence 
of 1.13 ± 0.05 x 10~ 6 erg cm -2 in the 15 — 150 keV band. The 1 second peak photon flux 
measured from T - 0.124 s in the 15 - 150 keV band is 12.1 ± 0.4 phones" 1 . All BAT 
errors are at the 90% confidence level. 

In the 15 — 350 keV lightcurve (Figure 1) at least 20 statistically significant peaks with 
FWHMs in the 5 — 15 ms range (Barthelmy et al. 2006) can be seen. There is no periodic 
structure in the lightcurve for at least the first 100 s. There is also no evidence of extended 
emission in the T + 1 to T + 300 s range at an upper limit of 0.001 counts detector -1 s -1 
(3<t). This is consistent with an upper limit on the flux ratio between the initial peak and 
the peak of any potential extended emission of ~ 2000. SHBs 050724 and 051227 had flux 
ratios of 46 and ~ 10, respectively (Barthelmy et al. 2006). 

A four channel lightcurve (Figure 2) reveals that the burst is a short hard burst. A 
lag analysis confirmed the burst to cleanly reside in the short hard burst class (Figure 3 in 
Norris & Bonnell 2006). The measured lags are 0.8 ± 0.6 ms [(50 - 100 keV)/(15 - 25 keV)] 
and 0.3 ± 0.7 ms [(100 - 350 keV)/(25 - 50 keV)]. 

2.2. KONUS-Wind Analysis 

KONUS-Wind (K-W) triggered on GRB 060313 at T (K - W) = 00 : 12 : 06.354 UT 
(Golenetskii et al. 2006). It was detected by the SI detector which observes the south 
ecliptic hemisphere; the incident angle was 57.9°. The propagation delay from Swift to 
Wind is 0.045 s for this GRB, therefore, correcting for this factor, one sees that the K-W 
trigger time corresponds to T — 0.085 s. 

The GRB time history was recorded in three energy ranges: Gl (21 — 83 keV), G2 
(83 — 360 keV), and G3 (360 — 1360 keV). Thanks to the high intensity of the burst, the 
K-W reserve system (so called 'nonius') was triggered providing 2 ms resolution record in 
the G2 and G3 ranges of the entire burst (up to T (K — W) + 0.762 s). Five spectra in 
101 channels were accumulated during the burst (the total number of spectra accumulated 
during the trigger record is 64). The light curve (Figure 1) shows several multi-peaked pulses 
with a total duration of approximately 0.8 s. As with the BAT multi-channel light curve, 
the three channel KONUS light curves (Figure 3) illustrate that the burst is a short hard 
burst. 
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2.3. XRT Analysis 

An uncatalogued, fading source was discovered during ground analysis (Pagani et al. 
2006b) of the XRT data. The X-ray afterglow was observed in Windowed Timing (WT) 
mode for 26 s. As the source faded the XRT switched into Photon Counting (PC) mode, thus 
providing 2-dimensional spatial information (for a complete description of the XRT modes 
see Hill et al. 2004). To avoid pile-up in the PC data during the first orbit of observations, the 
XRT lightcurve was extracted in the energy range 0.2 — 10.0 keV using an annulus centered 
on the afterglow position with an inner radius of 2 pixels and applying a PSF correction. The 
X-ray lightcurve (see Figure 4) during the first orbit of observation manifested variability, 
with spikes and fluctuations superimposed on the overall fading behavior. The initial decay 
slope was a = 1.25 ± 0.15. The lightcurve profile becomes smoother in the later orbits 
(> 4200 s after the BAT trigger) and the decay curve can be fitted by a power law with 
index of 1.46 ± 0.08. In addition, the flux at these later times is larger by a factor ~ 4 
compared with the extrapolation of the earlier X-ray lightcurve. The X-ray lightcurve is 
shown in Figure 4, with the count rate converted into 0.2 — 10.0 keV unabsorbed flux using 
the best spectral model fit. 

The X-ray spectrum from the Photon Counting data was extracted in the 0.3 — 10.0 keV 
energy range, with a binning of at least 20 counts per energy bin, using the latest versions 
of the ancillary and response matrix files. The PC spectrum of the first orbit of observation 
was extracted excluding the central 2 pixels to avoid pile up. The spectrum can be fitted 
by a simple absorbed power law, yielding a photon index of 1.53 ± 0.10 and an absorption 
column density Nh consistent with the Galactic value of 4.7 x 10 20 cm -2 (Dickey & Lockman 
1990). The spectrum softens during the later orbits: the best power law model fit of the 
afterglow spectrum extracted from 4.1 to 24 ks after the GRB trigger yields a photon index 
of 1.96 ± 0.09, fixing the N H at the Galactic value. 

2.4. UVOT Analysis 

The Swift spacecraft slewed promptly when the BAT detected GRB 060313, and UVOT 
began imaging the field shortly after the BAT trigger. The UVOT took 44 exposures of the 
field between 78 and 67 783 s after the BAT trigger. For the majority of these exposures the 
afterglow was too faint to be detected. The UVOT photometry that is used in this paper is 
presented in Table 1. The UVW2 detection (see below) implies that the Lyman limit must 
be blueward of ~ 250 nm, suggesting an upper limit on the redshift of z ^ 1.7. To further 
constrain the redshift, z, a power-law UV/optical SED with an unconstrained spectral slope, 
(3, as well as a fixed Galactic E(B — V) was assumed. The redshift was then varied and 
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the spectrum for a given (3 was modified according to the parameterization of the average 
Lyman absorption (Madau 1995). The magnitudes relative to each UVOT filter were then 
estimated for a grid of (3 and z, with a corresponding x 2 value for each (3-z pair. The grid 
of x 2 values were plotted to illustrate how /3 and z could be constrained. The probability 
distributions were also projected onto each axis to obtain the confidence intervals. Neither 
parameter is very well constrained, but the best (5 ~ 1.93 ± 0.22 and z < 1.1 at the 90% 
confidence level. The most likely redshift is z = 0.75, but it is minimally peaked. 

The UVOT has UBV filters which approximate the Johnson system and three ultraviolet 
filters: UVW1 with a central wavelength of A c = 251 nm, UVM2 with A c = 217 nm, and 
UVW2 with A c = 188 nm. Examination of the first settled observation (a 200 s exposure 
in V) revealed a new source relative to the Digital Sky Survey inside the XRT error circle. 
This source had a magnitude of V — 19.56 ± 0.32. The UVOT position of the source is 
RA = OAh 26m 28.429s and Dec = -10° 50' 39.13" (J2000), with an internal accuracy of 
±0"01 and an absolute astrometric uncertainty of / /56 (90% containment). This position is 
1"3 from the reported XRT position (Pagani et al. 2006b). It is inside the XRT error circle 
and consistent with reported ground-based detections (Levan & Hjorth 2006; Thone et al. 
2006). Subsequent exposures showed the source to be fading. 

We performed photometry on each UVOT exposure using a circular aperture with a 
radius of 2" centered on the position of the optical afterglow. This radius is approximately 
equal to the FWHM of the UVOT PSF. The PSF varies with filter and with spacecraft 
voltage, so we did not match the extraction aperture to the PSF for each exposure. The 
PSF FWHM, averaged over the voltage variations, ranges from l'/79 ± 0"05 for the V filter 
to 2 / /17±0"03 for the UVW2 filter. The background was measured in a sky annulus of inner 
radius 17" 5 and width 5" centered on the afterglow. 

Aperture corrections were computed for each exposure to convert the 2" photometry to 
the standard aperture radii used to define UVOT's photometric zero points (6" for UBV and 
12" for the UV filters). The aperture correction procedure also accommodates the variable 
PSF. Approximately six isolated stars were used to compute the aperture correction for 
each exposure. The RMS scatter in the mean aperture correction for a single exposure was 
typically ~ 0.07 mag. The RMS scatter for each exposure was added in quadrature to the 
statistical error in the 2" magnitude to obtain the total la error in each point. All detections 
above the 2a significance level are tabulated in Table 1. 

Since the UVOT is a photon-counting device, we have corrected all of our data for 
coincidence loss; however, the afterglow has V > 19, so coincidence losses are negligible, 
typically less than 0.01 mag. The zero points used to transform the instrumental UVOT 
magnitudes to Vega magnitudes are listed in Table 2. They are taken from the latest in- 
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Table 1. Swift/XJVOT photometry of the afterglow of GRB 060313 
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Table 2. Zero points applied to the photometry of the afterglow of GRB 060313 



Filter 


ZP 


Error 


V 


17.88 


0.09 


B 


19.16 


0.12 


U 


18.38 


0.23 


UVW1 


17.69 


0.02 


UVM2 


17.29 


0.23 


UVW2 


17.77 


0.02 


White 


19.78 


0.02 
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orbit measurements as obtained from the HEASARC Swift/UVOT Calibration Database 
(CalDB) 3 . Color terms were not applied to the photometric calibrations, but preliminary 
calibrations of on-orbit data suggest that they are negligible. 

The UVOT photometry was adjusted to the [/-band by assuming that the optical spec- 
trum was a power law with the same slope as the X-ray spectrum (f3 x = 0.96). The resultant 
UVOT photometry is provided in Table 1, Adjusted [/-band Mag column, and is shown in 
Figure 4. In addition, the UVOT photometry was also adjusted to the ?7-band assuming the 
cooling break is between the optical and X-ray bands at t > 1000 s. The assumed spectral 
slope in the optical is f3o = 0.46. Although the scatter looks slightly less for the f3o = 0.46 
case, the RMS residual for both cases is 0.45 mag (see below). 

As seen in Figure 4, there appears to be flaring in the UVOT lightcurve between 3000 — 
40 000 s. To determine the statistical significance of these potential flares, the shifted UVOT 
photometry was fitted with a single power law (Xred = 2-59 & rms sca tter = 0.45 mag for 
flo = 0.96; xied = 2-40 & rms sca tter = 0.45 mag for f3 = 0.46). Based on this fit, the 
fluctuations appear real. If the fluctuations represent flares and a single power law is fit to 
the bottom of the photometry data (i.e. a line is placed through the lightcurve while ignoring 
the three flares), the flaring becomes even more significant (Xred = 5.16 & rms scaUe r = 0.50 
mag for (3o = 0.96). A further check was made to determine if the fluctuations were caused 
by instrumental or other effects. Light curves for seven comparison stars (which were selected 
based on similar magnitudes to the afterglow's) were constructed in each filter after shifting 
to the [/-band as described above. For each comparison star a weighted mean magnitude 
and the residual about the weighted mean was computed. The RMS residual of all the 
observations of all the stars is 0.09 mag, which is much smaller than the RMS residual about 
the power law fit to the afterglow's light curve. We therefore conclude that the comparison 
stars do not show the same fluctuations as the optical afterglow and that the fluctuations 
are intrinsic to the afterglow. Although optical flaring has been reported for long bursts (cf. 
Holland et al. 2002; Jakobsson et al. 2004), this is the first short burst with optical flaring. 



2.5. Spectral Analysis 

Joint spectral analysis was carried out using the BAT data between 15 and 150 keV 
and the KONUS data from 20 to 3000 keV.The spectra were fit by a power law model 
with an exponential cut off: dN/dE ~ E~ r exp~( 2 ~ r ^ E / Ep where E p is the peak energy 
of the vE v spectrum and T is the photon index. A fit to the Band (GRBM) model was 



3 http: // swift .gsfc.nasa.gov / docs/heasarc / caldb /swift / 
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also attempted. No statistically significant high energy power-law tail was established in 
any fitted spectrum. The time-integrated spectrum is well fit with T = 0.61 ± 0.10 and 
E p = 947^3 keV (x 2 = 115/123dof). Using these parameters, the calculated burst fluence 
is 1.291q3i x 10" 5 erg cm" 2 and the peak flux is 5.99l5'59 x 10 -5 erg cm" 2 s" 1 , as measured 
from To = 0.475 s on a 16 ms timescale. All errors are at the 90% confidence level. Joint fits 
were also made to the same cut-off power law model for three specific time intervals denoted 
on Figure 3. The results are summarized in Table 3. Significant spectral evolution in both 
the low-energy photon index and E p is evident. 

We also constructed a spectral energy distribution (SED) of the optical and UV after- 
glow for 14 ks after the trigger. This time was chosen because the lightcurve is relatively 
well-sampled in the different UVOT bands around that point, requiring the minimum of 
interpolation in constructing the SED. For the B, U, and UVW1 bands respectively, a 
power-law was fitted to the data points between 4000 and 40 000 s and the count rates at 
14 ks were interpolated from the fitted curves. For the V, UVM2, and UVW2 bands, for 
which there was only one point per filter in the 4000 — 40 000 s time range, an average power- 
law decay fitted jointly to B, U, and UVW1 was renormalized to the single points in each of 
V, UVM2, and UVW2. The count rates at 14 ks in these filters were then determined from 
these extrapolated curves. The count rates from all six filters were used to make spectral 
files for fitting with XSPEC using the latest UVOT response matrices (version 102). The 
errors on the count rates consisted of a sum in quadrature of the photometric error and a 
10% systematic error from the calibration uncertainty on the filter effective areas. 

An XRT spectrum for the time interval 4100 — 24 000 s was extracted for joint fitting 
with the UVOT, and renormalized to the total count rate expected at 14 ks. Fitted by itself, 
the XRT spectrum was slightly softer than that obtained from the whole time interval, with 
a (Galactic absorbed) power-law slope of 1.97 ± 0.09. The XRT and UVOT spectra were 
then fitted simultaneously, with a model consisting of power-law absorbed by Galactic gas 
and dust as well as dust intrinsic to the GRB host. 

A good fit was obtained with the power-law with Galactic dust and gas alone; no dust 
intrinsic to the source was required. Since we cannot fit the redshift of any dust feature we 
are unable to use the SED fitting to constrain the redshift of the GRB in this case. The 
full parameters of the joint UVOT-XRT fit are given in Table 4, and the UVOT and XRT 
spectra are plotted with the model and the ratio of data to the model in Figure 5. 



Table 3. Parameters of a joint BAT-KONUS fit to a model consisting of a power law with 

an exponential cut off. a 
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1.18lg;i| 
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a Two sets of fit results are shown for each time interval (t). In the first 
row for each time interval, the joint fitting was done without renormalization 
between BAT and KONUS. In the second rows, the relative normalization 
(Irn) is allowed to be a free parameter. In the case in which is not sta- 
tistically consistent with unity, the difference is attributed to the differences 
in the low energy parts of the BAT and KONUS spectra. The photon index 
is T, E p is the peak energy, x 2 /dof is the goodness of fit, and Prob is the 
probability. 



Table 4. Parameters of a power-law fit to the UVOT-XRT spectrum with Galactic dust 
and gas plus dust absorption intrinsic to the GRB host. 



Component Parameter Values 

Power-law (7) 2.03 ± 0.04 

Normalization 5.5 ± 0.4 x 10~ 5 photons cm~ 2 s _1 keV -1 (@ 1 keV) 

Neutral Galactic gas absorption (Nh) 4.74 x 10 20 cm -2 (fixed) 

Dust type Milky Way 

E(B-V) 0.0625 mag (fixed) 

z (fixed) 

X 2 red 1-142 (34 dof) 
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3. Discussion 

3.1. Prompt Gamma- Ray Emission 

Although GRB 060313 is the hardest of the Swift bursts (based on the (50-100 keV)/(25- 
50 keV) ratio), it is just above the average hardness for BATSE short bursts (see Figure 6). 
Since the BAT lightcurve of GRB 060313 is bright and therefore the statistical significance 
is high, we can be confident that many of the multi-peaked pulses in the lightcurve are real. 
Most BATSE short bursts, even when as bright as GRB 060313, did not reveal as much 
structure as is evident in this burst. The burst fluence as measured in the 15 keV — 3 MeV 
band was 1.29^Q3i x 10~ 5 erg cm" 2 , providing an isotropic gamma-ray energy release of 
3.4 x 10 52 ergs for an assumed burst of redshift 1. 

The spectrum during the burst had a positive spectral index during the first 0.19 s of 
the burst (f u oc u im ) and the time averaged low energy spectrum was /„ oc //°- 29±0 - 07 . These 
spectra can be produced in the synchrotron-self-Compton model for gamma-ray generation 
(Kumar et al. 2006). The high degree of variability of the gamma-ray lightcurve is probably 
due to a highly variable outflow from the central explosion. 

3.2. Afterglow 

The X-ray lightcurve has a profile similar to the canonical profile described in the 
literature (Zhang et al. 2006) with the exception that segment I (the time immediately after 
the prompt emission) has a much shallower decay slope (~ 2x less than the canonical). This 
XRT decay slope of a = 1.25 ± 0.15 beginning ~ 100 s after the burst trigger is much less 
than the early decay slopes seen at similar times for some short bursts: GRBs 050724 (a = 7; 
Barthelmy et al. 2005b), 050813 (a = 2.05 ± 0.20; Racusin 2006), 051210 (a = 2.57 ±0.11; 
La Parola et al. 2006), and 051227 (a = 2.2 ± 0.2; Barbier et al. 2006). The decay slope is 
more consistent with the short bursts: GRBs 050509B (a = 1.20lo;{jl; Gehrels et al. 2005), 
050709 (a £ 1; Fox et al. 2005), and 051221A (a = 1.3 ± 0.7; Capalbi et al. 2005). 

The initial XRT lightcurve has flaring in the first 1000 s after the trigger. The spectrum 
during this period is harder (photon index of 1.53 ± 0.10) than the best spectral fit of the 
late afterglow (photon index of 1.96 ±0.09). Flares in the early lightcurve of the short GRBs 
050724 (Barthelmy et al. 2005b; Campana et al. 2006; Grupe et al. 2006), 051210 (La Parola 
et al. 2006), and 060121 (Levan et al. 2006) have also been seen. Flaring in GRBs has been 
attributed to long lasting activity by the central engine with internal shocks continuing for 
hundreds of seconds (Burrows et al. 2005b; Falcone et al. 2006a,b; Pagani et al. 2006c; Zhang 
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et al. 2006; De Pasquale et al. 2006). Although the decay slope of GRB 060313 is consistent 
with GRBs 050509B, 050709, and 051221A, no early flaring was seen in these three bursts 
(Fox et al. (2005) report the possibility of a flare at early times for GRB 050709, however, the 
fitting of the lightcurve at early times was difficult since only one data point was available in 
the first day after the burst). Between 1000 — 3000 s no data was obtained due to occultation 
by the Earth; therefore, it is unknown what the lightcurve profile is at this time. However, 
this region in the lightcurve is consistent with either segment II or V in Zhang et al. (2006), 
therefore, there is clear evidence that there is some kind of energy injection 

The optical and UV lightcurve is the earliest (< 80 s) optical/UV detection of a short 
GRB to date. The early UVOT decay slope, beginning ~ 100 s after the burst trigger, is 
a = 0.13 ± 0.28 which is much flatter than the XRT decay slope over the same period. 
In addition, unlike the XRT, the UVOT lightcurve contains no flaring during this period. 
At t > 3000 s, the UVOT lightcurve is steeper (a = 0.43 ± 0.13) and includes three flares 
superimposed on the overall lightcurve while the XRT lightcurve has no flaring. No other 
short burst has so far been observed to manifest optical flaring at such early times. 

The XRT spectrum during the first 1000 s of the burst was proportional to z/-°- 53±01 ° in 
the 0.3 — 10 keV energy band. The spectrum in this entire energy band was a single power- 
law meaning it did not contain any characteristic synchrotron frequency. During this same 
period, the XRT lightcurve fell off as t -125 signifying that the synchrotron peak frequency 
was below the XRT band or less than 0.3 keV. This indicates that the synchrotron cooling 
frequency was either less than 0.3 keV or above 10 keV. If the cooling frequency were to be 
below 0.3 keV then the XRT spectrum would be v~' p l 2 (where p is the electron power-law 
distribution index) giving a value for p = 1.06 ± 0.2 which is highly unlikely (it is much too 
small for the X-ray spectrum of GRB afterglows). We therefore conclude that the cooling 
frequency for t < 1000 s was above 10 keV (cf. Sari, Piran, & Narayan 1998; Chevalier & Li 
1999; Rol et al. 2005). An X-ray lightcurve due to shocks in a uniform density circum-stellar 
medium (CSM) is expected to fall off as t~ a8 , for a photon index of 1.53, instead of the 
observed t -1 - 25 decay. The observed lightcurve decay is consistent with a wind-like density 
stratification for the CSM. However, the subsequent decline of the X-ray lightcurve and the 
spectrum are seriously at odds with this possibility and therefore we do not consider it to 
be a viable solution. The steeper than expected fall off of the lightcurve could be due to 
some combination of a structured jet, radiative loss of energy, and decreasing microphysics 
parameters (e e & e^; cf. Panaitescu et al. 2006b). 

The flux in the 0.3 — 10 keV band at 100 s from the forward shock, assuming the XRT 
band below the synchrotron cooling frequency and a burst redshift of 1 which provides a 
luminosity distance (D 2 8) ~ 1, is (eq. B7 from Panaitescu & Kumar (2000) is used since we 
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assume Vi < v < v c and a homogeneous external medium) 

/, « 5 x 10~ 9 EH 4 n 1 / 2 ^ ^ erg cm" 2 s"\ (1) 

where E 52 is the isotropic equivalent energy in units of 10 52 ergs, e e -i = e e /0.1 (e e (p — — 
2) is the energy fraction in electrons), en-i is the energy fraction in the magnetic field, n is 
the density of the CSM, and we used p = 2.06 as suggested by the early X-ray spectrum 4 , 
in calculating the numerical values in the above equation. Note that the error in the burst 
redshift has a very small effect on the X-ray flux formula given above as f x has almost a 
linear dependence on the isotropic burst energy and therefore the effect of uncertainty in 
distance on _E 52 and f x nearly cancel each other. Substituting the observed value of the 
X-ray flux of 3 x 1CT 10 erg cm -2 s _1 in the above formula we find 

^2 /4 n 1/2 6 e ,- ie | /4 _ 2 -6xl0- 2 . (2) 

Moreover, applying the constraint that the cooling frequency is greater than 10 keV at 10 3 s 
(because the photon spectral index is 1.53) we obtain 

^ 5 1 2 /2 ^4 /2 -2<5xl0- 3 . (3) 

Equation (3) follows from the expression of cooling frequency given in the literature (eq. 27 
from Panaitescu & Kumar (2000) assuming K <C 1). Combining equations (2) and (3) we 
find: e e _i E 52 > 1- Moreover, for p = 2.06 the fraction of energy of the shocked fluid in 
electrons is equal to 1.76e ej _i which cannot be greater than 1, and therefore e e _i < 0.6 and 
E 52 > 1-8. The 15 keV — 3 MeV 7-ray fluence suggests that the isotropic equivalent energy 
in photons, assuming z — 1, was E52 = 3.4, which is consistent with the above derived 
constraint on the energy in the explosion of E52 > 1-8. Since the radiative efficiency for the 
prompt 7-ray emission must be less than 100%, and likely of order a few tens of percent, 
this suggests E 52 ~ 6, and e 6) _i ~ 0.5. Substituting these values in equation (2) we find 

™ 1/2 4 /4 -2~lxl(n 2 - (4) 

We next make use of constraints on the synchrotron injection frequency (the synchrotron 
frequency at which the bulk of the electrons which are injected by the shock front radiate) 
which is given by (e.g. eq. 22 in Panaitescu & Kumar 2000) 

*i = 306 Eg ef_ 2 1~ 3/2 eV. (5) 



4 Since /„ oc v 53 , then p — 2.06 ± 0.2. The lower bound (1.86) would lead to a negative energy fraction 
of electrons which is an unphysical event. 
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Since the XRT lightcurve is falling off as t~ L25 starting from the very first XRT data point 
at t 2 — 1 (where t 2 = t/100 s), we conclude that z/j < 0.3 keV at this time. If z/j > 0.3 keV, 
then the XRT lightcurve would be rising with time as t l l 2 (cf. Panaitescu & Kumar 2000). 
Using equation (5) for z/j we therefore conclude that 63,-2 < 2.6; E 52 = 6 & e e> _i = 0.5 were 
used in this estimation. Moreover, at t = 1 ks or t 2 = 10, vi should be greater than 2 eV 
since the UVOT lightcurve is flat during the time interval 100 — 1000 s. This constraint, 
and equation (5), leads to €3-2 > 0.1. From these upper and lower limits on €3-2 it seems 
reasonable to infer that €3-2 ~ 1) which when substituted into equation (4) gives the density 
of the CSM to be n ~ 10~ 4 cm" 3 . This is consistent with the combined UVOT/XRT SED 
which reveals that no intrinsic dust is required to model the environment of the burst. This 
appears that the burst went off in a low density medium that is usually found in bubbles or 
the outskirts of galaxies. 

If the UVOT band lies below u i: the lightcurve is expected to rise as t 1 ' 2 (cf. Piran 1999; 
Panaitescu & Kumar 2000). However, the lightcurve seems to be nearly flat for the initial 
10 minutes or so, thus, we conclude that z/j is below the UVOT band. It is likely that the 
same process that caused a steeper fall off of the early X-ray lightcurve, by ~ t ' 45 , could 
have also flattened the rise of the early UVOT lightcurve from the expected t 1 ^ 2 to The 
general decline of the optical lightcurve between 10 and 40 ks is roughly the same as the 
early X-ray lightcurve and that is consistent with the forward external shock model. 

The X-ray lightcurve for t > 4 ks is steeper than the early lightcurve by t 1 / 4 and that 
is almost certainly a result of the cooling frequency dropping below the XRT band; the 
spectrum during this stage also became softer as expected for a cooling transition. However, 
the flux at t > 4 ks is larger by a factor ~ 4 compared with the extrapolation of the earlier 
X-ray lightcurve. This suggests quite a substantial amount of energy being added to the 
forward shock, almost twice the amount of the initial energy in the explosion, during the 
time period of 1 to 4 ks. It is also possible that the enhanced X-ray flux could be due to the 
GRB ejecta running into a denser shell of CSM. 

There are fluctuations in the UVOT lightcurve between 4 and 40 ks whereas during 
this period the X-ray lightcurve has a smooth powerlaw decline. This interesting behavior 
is not expected for energy being added to the decelerating forward shock (due to late time 
central engine activity, for instance) whereas this is consistent with fluctuations in the den- 
sity of the CSM. The observed flux at frequencies above the cooling frequnecy (z/ c ) is very 
insensitive to the CSM density structure and this could explain the smoothly declining XRT 
lightcurve. The flux at frequencies between z/j and u c varies with density as n}l 2 (cf. Piran 
1999; Panaitescu & Kumar 2000), and hence density fluctuations will be reflected in the 
optical lightcurve if u opt < v c \ the amplitude of the UVOT flux variation corresponds to a 
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factor of about 2 variations in CSM density, and the short timescale for the variability sug- 
gests small scale fluctuation in the CSM. We note that the early X-ray lightcurve (t < 10 3 s), 
when the XRT band was below u c , also showed fluctuations, lending additional credence to 
the interpretation of CSM density fluctuation. 

The variations in the XRT and UVOT lightcurve occur on very short timescales (dt/t <C 
f, where dt is the variability timescale). These small amplitude fluctuations can arise in the 
external shock with dt/t < 1; for instance, when the shock front crosses a density clump of 
size much less than r /Y (r is the radius of the external shock and Y is the Lorentz factor) it 
will increase the observed flux a bit. This increase will last for a time dt < t where dt/t is 
determined by the clump size, and the amplitude of fluctuation by its mass. 

This short burst offers conflicting evidence as to the nature of the progenitor and its 
surrounding medium. On the one hand the afterglow data suggests that the burst went off 
in a low density medium and on the other hand fluctuations in the UVOT lightcurve and a 
lack of corresponding fluctuation in the X-ray data suggests that the CSM density varied by 
a factor of order 2 on a length scale of a few times 10 17 cm. 

The low density of the circum-stellar medium has fluctuations that one normally as- 
sociates with the surroundings of a massive star and inhomogeneities of a windy medium. 
However, the prospect of the progenitor being a massive star is unlikely. The mean density 
of the medium for this burst is small (~ 1CT 4 particles cm -3 ) for a massive star. In addition, 
although it can not be ruled out, no host galaxy was found at the location of GRB 060313. 
Instead, a search in the extragalactic databases indicates that there are 12 galaxies within 
5' of GRB 060313, of which 6 are within 2' and one within 30" (B042408.28-105721.7). Could 
this possibly be a galaxy cluster with which the GRB is associated? Currently there are no 
known redshifts or color information on these galaxies; therefore, it can not be determined 
for certain. However, if we assume that there is no host galaxy at the position of the GRB 
and that it is part of this potential galaxy cluster, it is unlikely that a massive star can 
travel such great distances. Massive stars live for about 10 million years, and their peculiar 
velocity is not more than a few tens of kms" 1 . Therefore, in their lifetime, a massive star 
can only travel a distance of less than a few hundred pc. 

In contrast, the natal kick of a neutron star (~ 200 — 1000 kms -1 ) and the typical 
life times of compact binaries (~ 0.1 — 1 Gyr) allow for DNS and NS-BH binaries to travel 
great distances. Assuming a redshift for GRB 060313 of z = 0.16 (for GRB 050709 - the 
closest measured redshift for a short GRB; Fox et al. 2005), z = 0.35 (the average measured 
redshift for short GRBs; La Parola et al. 2006), z = 0.55 (for GRB 051221A - the farthest 
measured redshift for a short GRB; Soderberg et al. 2006), and z = 1.1 (for the UVOT upper 
redshift limit), a compact binary progenitor could have been associated with a galaxy that is 
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^ 6.2', 3.5', 2.7', and 2.1' away, respectively. These distances are consistent with the distances 
to some of the galaxies in the potential cluster. Without the early lightcurves, particularly 
the UV/optical lightcurve, we would have automatically assumed that the progenitor of this 
GRB was a compact binary merger. However, the data suggests that GRB 060313 may arise 
from a mechanism other than the established compact binary model. 

An alternative possibility is that the rapid variabilities in the UVOT lightcurve are 
produced by late-time central engine activity similar to those produced in late-time X-ray 
flares seen in many other GRBs (e.g. Burrows et al. 2005b; Zhang et al. 2006; Levan et al. 
2006). Such a possibility has been raised to discuss other rapid optical variabilities of pre- 
Swift afterglows with At < t (e.g. Ioka, Kobayashi, & Zhang 2005). In this scenario, no 
CBM density fluctuation is required. The lack of variability in the X-ray band may be due 
to the fact that the EpS of these late flares are very low (i.e. in the UV/optical regime), 
so that their contributions to the X-ray band are negligible compared with the afterglow 
level. Within the late internal shock scenario, soft flares are in principle possible given the 
combination of a low luminosity and a large Lorentz factor as a result of less baryon loading 
at later times (Zhang & Meszaros 2002). The challenge is how to restart the central engine 
at very late epochs in short GRBs. The same problem has also been encountered when 
interpreting the multiple late time flares of GRB 050724 (Barthelmy et al. 2005b), and some 
suggestions have been made to explain its behavior (e.g. Barthelmy et al. 2005b; Perna, 
Armitage, & Zhang 2006; Proga & Zhang 2006; Dai et al. 2006; Faber et al. 2006). 



4. Conclusion 

The multi-wavelength lightcurve of the SHB GRB 060313 reveals a multi-peaked struc- 
ture in the 7-ray, hard X-ray, X-ray, and UV/Optical bands at different epochs of the prompt 
and afterglow emissions. The hard X-ray spectrum during the prompt emission had a posi- 
tive spectral index consistent with the synchrotron-self- Compton model for 7-ray generation. 
The large number of peaks in the 7-ray lightcurve is possibly the result of the central engine 
generating an extremely variable outflow. 

The early (70 — 1000 s) X-ray temporal decay slope (a = 1.25 ± 0.15) is exceptionally 
steep for the observed X-ray spectrum and manifests flaring components superimposed on 
the lightcurve. During the same period the UV/Optical lightcurve reveals a flat decay profile 
(a = 0.13 ± 0.28) with no flaring. This behavior can be attributed to a combination of a 
structured jet, radiative loss of energy, and decreasing microphysics parameters. The X-ray 
spectrum indicates that the cooling frequency was not much greater than 10 keV at 1000 s 
after the burst. At later times (> 3000 s) the UV/Optical lightcurve manifests flaring on top 
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of a decaying profile (a — 0.43 ± 0.13) while the X-ray exhibits no evidence of flaring on top 
of a lightcurve that is steeper than at earlier times. An interpretation of the data points to 
a burst that occurred in a low density medium and that the CSM density varied by a factor 
of order 2 on a length scale of a few times 10 17 cm. Although the collapse of a massive star 
as the progenitor mechanism is doubtful (though it can not be ruled out), the favored model 
of a compact binary merger also creates problems in explaining the burst's behavior. One 
explanation for this behavior is that at least some short bursts may be the result of a different 
mechanism other than the traditional compact binary model. An alternate explanation is 
that late-time central engine activity is injecting energy into the UV/optical regime but not 
into the X-ray regime. 

This is the first short burst that has manifested this kind of behavior. Because the 
sample size is small, more observations of these objects are needed in order to determine 
whether GRB 060313 is the norm or anomalous for short GRBs. In addition, more detailed 
theoretical modeling is needed to establish a clearer picture of the mechanism. 
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Fig. 1. — BAT and KONUS-Wind lightcurves. The BAT 15 — 350 keV lightcurve (upper 
panel) has at least 20 statistically significant peaks with FWHMs in the 5 — 15 ms range. 
The average BAT error is ^SOcountss" 1 . The KONUS-Wind 21-1360 keV lightcurve (lower 
panel) also exhibits several multi-peaked pulses with a total duration of approximately 0.8 s. 
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Fig. 2.— BAT four channel (15 - 25, 25 - 50, 50 - 100, & 100 - 350 keV) light curves. The 
average errors are 655, 670, 600, and 530 countss -1 , respectively. The fifth panel from the 
top shows the photon power-law index for a simple power law fit. The bottom panel shows 
the fluence hardness ratio S3/S2 between the 50 — 100 keV and 25 — 50 keV bands, where the 
fluence was calculated based on a power-law fit. The hardness ratio for the time integrated 
spectral fit is 2.43. 
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Fig. 3. — KONUS-Wind lightcurve of GRB 060313 in three energy bands. The two bottom 
panels show the hardness ratios G2/G1 and G3/G2 which demonstrates significant spectral 
evolution during the burst. The dashed vertical lines denote the boundaries of the three 
intervals over which we extracted the spectra in Section 2.4. 
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Fig. 4. — Combined BAT, XRT, and UVOT lightcurves. The green represents the BAT 
values extrapolated into the XRT energy range, the red represents the XRT 0.3 — 10 keV 
fluxes, and the blue represent the UVOT values adjusted to the [/-band. The two insets 
are provided to show the flaring of the UVOT (insert (a)) and X-ray (insert (b)) lightcurves 
after 3000 s and during the first 1000 s, respectively. 
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Fig. 5. — Top panel: Combined UVOT and XRT spectral energy distribution (SED) of 
the afterglow at 14 ks after the burst (red = V, green = B, dark blue = U, light blue = 
UVW1, pink = UVM2, yellow = UVW2, black = X-ray). The model fit is over plotted; 
this model, which was fitted to both the UVOT and XRT data simultaneously, consists of 
a single power-law with a slope of 2.03 ± 0.04 and Galactic neutral gas absorption and dust 
extinction (the full model parameters are listed in Table 4). This power-law is softer that 
that fitted to the complete XRT dataset (which is biased towards earlier times due to the 
declining count rate). No absorption or extinction intrinsic to the GRB host is required in 
the fit. Bottom panel: Ratio of this model to the UVOT and XRT data. 
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Fig. 6. — Hardness ratio versus T 90 for the BAT and BATSE bursts. This plot shows how 
the spectral hardness and duration of GRB 060313 compares to previous BAT bursts (X 
symbols), and BATSE bursts (dots). One can see that GRB 060313 is the hardest burst yet 
detected by BAT, although it falls well within the BATSE hardness distribution. 



